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Abstract . — Galactic open clusters provide an abundant sample of stellar aggregates of various sizes, ages and metal 
abundances, apt to constitute a template for comparison with star systems in other galaxies. In this paper we present 
and discuss a standard methodology to synthesize U,-B,V fluxes and colours, and apply it to a set of 138 open clusters. 
Results are compared with previous available ones in the literature. We were able to calibrate a mass-luminosity 
relation by which we evaluated the mass of ~ 400 open clusters, leading to a well deflned present-day mass function. 
The number-complete sample of galactic open clusters presented in Battinelli & Capuzzo-Dolcetta (1991) is enlarged 
of a 15%. 
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1. Introduction 

It is well known that in our Galaxy the two families 
of open and globular clusters are quite well distinct in 
sizes, spatial distribution, metallicity content, kinemati- 
cal properties and even ages. In the galaxies of the Local 
Group, such a sharp distinction is no longer possible: as 
an example, the popoulous family of stellar clusters in 
the Magellanic Clouds spans with continuity a huge inter- 
val of masses, ages, heavy elements content, without any 
evident distinction in classes. Therefore, any meaningful 
comparison between stellar systems in galaxies requires a 
careful definition of the sample characteristics. 
Moreover, the data relative to crowded star clusters in 
galaxies are mainly relying on their integrated light. As a 
consequence, the possibility of a valid comparison between 
stellar populations in star clusters of different galaxies 
stands also on an uniform and objective way to com- 
pute integrated fluxes and colours for a reliable sample 
of galactic open clusters. In this regard, we note that, in 
the literature, various sets of integrated fluxes and colours 
for galactic open clusters have been presented (see, e.g.. 
Gray 1965; Sagar, Joshi & Sinvhal 1983, hereafter SJS; 
Spassova & Baev 1986, hereafter SB; Lynga 1987; Pandey 
et al. 1989). Of course, it is important to know precisely 
the way followed to compute the integrated fluxes; un- 
fortunately the richest sample of data (Lynga 1987) is 
not accompanied by any description of the method. In 
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particular, it is quite clear the crucial importance of the 
selection of bright stars as member of a cluster when de- 
termining its integrated magnitude. 

Once the relevant integral quantities are computed, an- 
other crucial point (already addressed in Battinelli & 
Capuzzo-Dolcetta 1991, hereafter BCD) is the statisti- 
cal completeness of the sample of clusters to use as a 
comparison template. 

The main aim of this paper is to present and dis- 
cuss a standard method to reliably evaluate integrated 
magnitudes and colours (and some derived quantities) for 
galactic open clusters; the FORTRAN package able to 
produce synthetic U,B,V, fluxes following the method de- 
scribed in Sect. 2.1. This package can be easily adapted to 
deal with input data coming from different sources; this 
is of particular importance due to the rapidly increasing 
availability of digitized data in astronomical data centers 
(see for instance the open cluster data by Mermilliod 1986, 
available at the Strasbourg data Centre) 



2. Integrated fluxes 

2.1. The standard method of computation 

Our methodology to evaluate U,B, V integrated fluxes for 
galactic open clusters is simple; using the most suitable 
and up-to-date data (HR diagrams, radial velocities and 
proper motions when available) for open clusters available 
in the literature, we try to define a standard procedure for 
membership determination and for flux summation. Such 
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a procedure (described below) has been automatized, after 
the necessary tests. 

Standard assumptions made throughout this work are: 
the reddening is the same for all the stars of a given clus- 
ter, and was taken from Lynga 1987; the ratio Av/Eb-v 
is set to 3.1; when Ejj-b is not available, the relation 
Ejj^B = 0.72Eb-v is assumed; when the data available 
refer to the RGU system (Becker 1946) they are trans- 
formed to U BV via the relations given by Steinlin (1968) 
(this happened for only 10 clusters). 

Regarding to membership, we relied on quoted refer- 
ences and decided to use criteria: i) to exclude stars just 
when explicitly suggested by radial velocity or proper mo- 
tion in the source paper; ii) to exclude stars external to 
a certain angular field, as given in the reference quoted, 
which is usually defined on the basis of stellar counts; 
iii) to exclude very bright blue stars just when red gi- 
ant are surely members, and just when they are at least 
3 mag brighter than the turn-off (this to avoid to exclude 
possible blue-stragglers, which have been found in many 
clusters). Anyway, it cannot be ruled out the possibility 
that a population of blue bright stars is a real feature 
of the cluster, even when the presence of evolved red gi- 
ants seems in contradiction to that. Indeed, for example, 
bright blue stars can be part of a second generation. For 
this reason, the exclusion of blue bright stars is done only 
when they arc very few (< 2). 

In addition to the membership determination, the 
main sources of error in the evaluation of the absolute 
synthetic magnitudes arc: the photometry of the single 
stars (standard errors of ~ 0.02 and 0.06 mag, for photo- 
electric and photographic data, respectively); the distance 
modulus determination, whose typical error is < 0.4 mag 
(see SJS); the reddening, carrying an indeterminacy of 
±0.1 mag (see SJS), when the hypothesis of constant 
Eb-v inside the cluster is acceptable; the incoming in- 
completeness of the HR diagram of the cluster at faint 
magnitudes, whose importance is weighted by the slope 
of the luminosity function. 

Globally, we can confirm the SJS conclusion that the 
maximum error in magnitude is about 0.5 m,ag and about 
0.2 mag in colour. 

2.2. The comparison of our results with 
available data 

The Lynga's Catalogue of Open Cluster Data (1987) con- 
tains values of integrated V and B — V ior k, 400 open 
clusters, mostly computed by Skiff. Unfortunately this 
data-base of synthetic fluxes, which is the most abundant 
available, has not been accompanied by a description of 
the methodology used and sometimes it is even difficult 
to identify the original source of data. Anyway, the abun- 
dance of Lynga's data suggests to check them with fluxes 
and colours obtained with an independent method whose 



reliability can be evaluated. In this regard, we decided 
to perform the comparison on a representative sub-set of 
the Lynga's sample consisting of clusters spanning large 
intervals of age, distance, and richness. 



Figure 1. Comparison between our (abscissa) and Lynga's 
integrated My (panel a) and {B — V)q (panel b). The sample 
of 53 clusters is composed only by clusters whose data are 
taken from the same reference source used by Lynga. Squares 
indicate clusters of the BCD complete sample. 

Figure 1 shows the comparison between our and 

Lynga's integrated magnitudes and colours for a sam- 
ple of = 53 clusters adopting the same sources of 
data. The agreement is good (the slope of the My cor- 
relation is a = 0.99 ± 0.03, with a correlation coefficient 
r = 0.98; a = 0.92 ± 0.05 and r = 0.93 for {B - V)o), par- 
ticularly for the bright and nearby clusters of the BCD 
sample (My < —4.5 and projected heliocentric distance 
d < 2 Kpc). In both our and Lynga's determination some 
(anomalous) too blue clusters are present. The reason for 
that is both an excessive average value of E{B — V) as 
deduced from Lynga's catalogue and the dispersion of 
individual reddenings around the average. Actually, we 
checked that the integral (B — V)o of Boc 7 (the only 
cluster of the 7 clusters bluer than {B — V)o = —0.5 
for which these data are available) changes from —0.6 to 
—0.24 when using individual E{B — V). This is due to a 
systematically small (respect to the average) reddening of 
the luminous stars in Boc 7. 

A systematic overestimate of our integral V-luminosity 
with respect to Lynga is clear at faint magnitudes, and it 
has been identified as due to the exclusion by Skiff of some 
very bright stars. These exclusions do not seem to be mo- 
tivated by a unique criterion for all clusters. Obviously, 
the exclusion of bright stars is more important for the 
computation of the integral magnitude of poor clusters. 
The difference in the integrated colours is not systematic, 
depending on the colour of the excluded stars. 

The inclusion of more updated individual cluster data, 
with respect to those used by Skiff, obviously increases 
the scatter in the integrated fluxes correlation as shown 
in Figure 2. 
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Figure 2. As in Fig. 1 for clusters of Table 1. 

A comparison of our integrated V magnitude, B — V 
and U — B colours with those given by Gray (1965), SJS, 
SB and Pandey et al. (1989) is given in Figure 3 and 4; 
the agreement is not as good as in the case of Lynga's 
data, due to both different references and methodology. 



Table 1 shows our integrated U,B,V magnitudes and 
colours for a sample of clusters composed by all the objects 
in Lynga's catalogue within the projected distance (onto 
the galactic plane) d — 2 Kpc from the Sun ( thus to 
enlarge the BCD sample) and of a selected set of clusters 
apt to check our method versus Lynga's data. 
The HR diagrams of the various clusters, as available 
from the literature, arc obviously not homogeneous in 
their completeness characteristics. For this reason, we have 
introduced Vum as the faintest magnitude of stars included 
in our flux summation. Of course, the importance of Vum 
is determined by the luminosity function of the cluster, 
in the sense that a variation of Vum is more relevant for 
clusters poor in high luminosity stars and/or having a 
steep slope of the luminosity function in the low main 
sequence. An estimate good within < 0.05 mag in the 
integral V magnitude (corresponding to ALy/Lv — 5%) 
requires to include in the computation of the flux, stars 
in an apparent magnitude interval at least AVmin = 6.0 
starting from the brightest member. 

2.3. The new complete sample 

The data base of all clusters with known distance and 
integrated My allows us to update the number-complete 
sample described in BCD. If we make the hypothesis 
of uniform averaged distribution of clusters (i.e. the as- 
sumption of a constant number density of clusters within 
cylinders of radius d centered to the Sun, see BCD Sect. 
2) 



Figure 3. Comparison between our integrated Aly (abscissa) 
and those computed by Gray (1965) (squares), Sagar et al. 
(1983) and Pandey et al. (1989) (crosses), and Spassova & Baev 
(1985) (triangles). 



Figure 5. Count of clusters brighter than My = —4.5 inside 
cylinders of radius d (in Kpc) centered at the Sun. The dashed 
straight line has slope 2, to correspond to a uniform distri- 
bution, and is normalized to the observed number of clusters 
within d = 1 Kpc. 



Figure 4. {B — V)o and {U — B)o correlations (panel a and 
b) of ours and Gray (1965), Sagar et al. (1983), Pandey et al. 
(1989) and Spassova & Baev (1985). Symbols as in Figure 3. 



Figure 5 shows how the sample of N = 115 clusters 
brighter than My = —4.5 and with projected distance 
d < 2 Kpc is number-complete. Hereafter we will re- 
fer to it as the 'complete' sample, whose relevant char- 
acteristics are given in Table 2. The present-day clus- 
ter formation rate as deduced by the counts of clus- 
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ters younger than Log t{yr) = 6.5 in the new sample 
is i? = 48 X 10^^ K pc^'^yr^^ , i.e. unchanged (within the 
10% error estimated in BCD) with respect to the BCD 
evaluation. 



2.4. The hyades cluster 

The Hyades are known to constitute a reference cluster 
apt to give a reliable fiducial unreddened main sequence 
and two-colour diagram. 

By mean of the HR diagram (187 members) given by 
Hagen (1970) wc get V = 0.49, My = -2.91, B - V = 
{B — V)q = 0.40; these data are almost coincident with 
Lynga's values {V = 0.50, My = -2.91, [B - V)o = 
0.40) in spite of the different number (380) of member 
stars. This means that ~ 50% of stars in Lynga's sample 
contribute negligibly to the total light. 
Adopting more recent data (Pels, Oort & Pels-Kluyver 
1975; Griffin et al. 1988; Gunn et al. 1988) we collected 256 
stars as certain members (on the basis of radial velocities), 
and obtained (with (M — m)v ~ 3.24 mag from Gunn et 
al. 1988) V = 0.56, My = -2.68, {B - V)o = 0.40. 



Figure 6. Colour-magnitude diagram for the whole set of ours 
and Lynga's clusters. The lower and upper solid lines are the- 
oretical evolutionary tracks (see text) of a solar composition 
cluster of total mass 25 Mq and 40, 000 M©, respectively. The 
thick curve is the location of the luminosity class V. 



3. Photometric results 

The results of Sect. 2.2 indicate the good degree of relia- 
bility of Lynga's integrated magnitudes. For this reason, 
the global results obtained from his B,V integrated pho- 
tometry stand almost unaltered by our enlargement of his 
sample of ~ 10%. Anyway, our computation of U fluxes 
for the in Table 1 suggested us to examine the photometric 
characteristics of our sample. 

Figure 6 is the integrated colour-magnitude diagram 
of the clusters of Lynga's catalogue and of Table 1 sam- 
ple. This figure seems to confirm the reliability of cluster 
integrated photometry because just few clusters lie below 
the luminosity class V (taken from Allen 1973). 

The colour-colour diagram in Figure 7 shows a defined 
trend with a slope in good agreement with the Z = 0.02 
theoretical model taken from Battinelli and Capuzzo- 
Dolcetta (1989). 

Note that in Figures 6 and 7 some clusters are present with 
anomalously blue {{B — V)o < 0.5) colour. Five clusters 
out of six arc poor and distant; the reason for their too 
blue colours is both an excessive average value of E{B — 
V) as deduced from Lynga's catalogue (or, equivalently, 
the photometry has a blue bias) and the dispersion of 
individual reddenings around the average. Actually, we 
checked that the integral (B — V)o of Boc 7 (the only 
cluster of the 6 for which these data arc available) changes 
from —0.6 to —0.24 when using individual E{B — V). This 
is due to a systematically small (respect to the average) 
reddening of the luminous stars in Boc 7. 



Figure 7. This colour colour diagram refers to clusters of Ta- 
ble 1. The thin curve is the theoretical model, while the thick 
one refers to luminosity class V. 

Figure 8 shows the My vs age relation; let us note that 
clusters lie between the theoretical curves corresponding 
to a 25 Mq and 4 x 10'' Mq cluster, as obtained using the 
theoretical model above. 

As it happens for Magellanic Cloud clusters (see Bat- 
tinelli & Capuzzo-Dolcetta 1989) we find that {U — B)q is 
slightly better correlated to the cluster age than {B — V)o, 
the least-square relation being: 

{U - B)o = OAGSLog t - 3.95, (1) 

the correlation coefficient is r = 0.643. 

4. A photometric estimate of open cluster masses 

The best way to evaluate the total mass of a stellar 
system is, in principle, given by dynamical considerations 
based on velocity and position data of member stars. 
However, suitable radial velocity and proper motion data 
are available just for 9 clusters, and so it may be helpful 
to resort to photometric data. A photometric estimate of 
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found |Am/m| ~ |AM/M|<10% when Y varies and Z 
is kept constant, and \ Am/m\ ~ |AM/M| <20% when Z 
varies and Y is kept constant. In some cases it was needed 
to extrapolate the m — L relation outside the 0.7 — 9Mq 
interval. A high mass extrapolation was required for 4 
clusters only; however, because of the scarcity of these 
massive (m > 9Mq) stars, the evaluation of the total 
cluster mass is not sensibly affected. An extrapolation 
below 0.7 Mq was necessary just for Hyadcs, where inde- 
pendent estimates of the mass allow to check the error. 
The bolometric corrections necessary to convert visual 
into bolometric luminosities have been taken from John- 
son (1966) and Morton and Adams (1968). Note that the 
bolometric correction induces an error small in compari- 
son with that due to other parameters (e.g. the distance 
modulus). With regard to the influence of a variation of 
the low-luminosity cut-off, Lvum, we checked that to 
have a relative error less than 5% in the cluster mass es- 
timate, all the stars in the range of at least 9 mag fainter 
than the brightest cluster star should be included. 



Figure 8. Panel a: integrated My vs age relation for the same 
sample of Fig. 6. Panel b: integrated {U — B)o vs age for the 
same sample of Fig. 7. The curves are the same theoretical 
models of Fig. 6. 



the total mass of a cluster is (see also Pandey, Bhatt & 
Mahra 1987, hereafter PBM) 



M = mil){m)dm c± / mijj{m)dm = 
= / m{Ly)(j){Lv)dLv 



(2) 



(where '^(m) and <p{Lv) are, respectively, the mass and 
^-luminosity functions of the cluster, and is the 

mass corresponding to Lviim)- Of course, this evaluation 
of cluster mass stands on a suitable mass-luminosity rela- 
tion: moreover, Lvurn should be low enough to allow the 
convergence of the last integral in Eq. (2) (hereafter mass 
and luminosities are in solar units). 

The approximations adopted are: i) to main sequence 
stars we assigned their ZAMS mass; ii) to post-main se- 
quence stars we assigned the turn-off mass value. The 
relative error introduced is evaluated to be less than 
10%. The mass-luminosity relation on the ZAMS has 
been calibrated on the basis of Mengel et al. 1979 
(0.7 < tti/Mq < 3) and Becker 1981 (3 < m/Mg < 9) sets 
of data suitably matched in the region overlap. A solar 
composition {X = 0.70, Y = 0.28, Z = 0.02) was chosen; 
this required interpolations in Y between 0.2 and 0.3 and 
in Z between 0.01 and 0.04 of Mengel et al. data. To give 
a rough evaluation of the sensitivity of the mass on the 
ZAMS at varying Y and Z in the mentioned intervals we 



Figure 9. Correlation of our mass evaluations (abscissa) and 
those by Schmidt (1963) (small black squares), Pandey et al. 
(1987) (black squares), and Bruch & Sanders (1983) (empty 
squares). Masses are in solar units. 



Figure 9 shows a comparison between our mass evalu- 
ations and those by Schmidt (1963), PBM, and Bruch & 
Sanders (1983). Our values are systematically lower than 
Schmidt's ones (which, being dynamical estimates, take 
also into account non luminous mass) and, in the average, 
slightly larger than PBM ones. Probably, PBM values are 
smaller than ours due to incoming incompleteness at their 
fainter magnitudes causing an apparent convergence of 
the total mass as given by Eq. (2). 

The slope of the least square fit to the relation between 
Schmidt's and our masses is close to 1: this indicates 
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(even with a large dispersion) the remarkable result of a 
roughly constant fraction of non-luminous mass to the 
total in open clusters. 

For the sake of comparison, Figure 10 shows our and 
PBM cumulative masses in func;tion of the individual mass 
for the Hyades cluster, for which we used the detailed set 
of data by Griffin et al. (1988). PBM give 7SMq for the 
total luminous Hyades mass; we obtain 256±3A'/0, against 
the 692Mq given by Schmidt and 390Mq by Gunn et al. 
1988, who both consider non-luminous matter. 



Figure 11. Logarithmic mass-luminosity relation for clusters 
of Table 1. The least square fitting relation is shown. Masses 
are in solar units. 



Figure 10. Value of the mass computed via Eq.(2) in func- 
tion of for the Hyades cluster. Crosses: this paper; dots: 
Pandey et al. (1987). Masses are in solar units. 



Figure 10 is an example of how the existence of a 
plateau in the cumulative mass frequency may be due to 
incompleteness, rather than to the actual convergence of 
the integral in Eq (2). 

The mass-luminosity relation for clusters in our sample is 
shown in Figure 11; the least-square fit to this relation is: 



LogM = -0.15Mv + 1.43, 



(3) 



with correlation coefficient r = 0.69. The standard devia- 
tion of data from the estimate is cr = 0.29, corresponding 
to a factor « 2 in the mass. 

An increasing trend of the M/Ly ratio with age is 
found even with a large dispersion (see Figure 12); the 
corresponding least square fit 



Log^ = 0A39Logt - 5.09, 



(4) 



with r=0.67. 

Individual values of cluster mass deduced via Eq. (4) 
are suitable just at the level allowed by the large dispersion 
in the data; anyway, it can be reliably used to obtain the 
statistical information given by the global mass function. 
The resulting mass histogram of all clusters in Lynga's 
catalogue with available age and luminosity is shown in 
Figure 13. 

The exponent of the least-square fitting to the decreas- 
ing part of the mass function is — 2.04±0.11 (r = 0.96) and 



Figure 12. M/Ly ratio vs age (logarithmic scales) for our 
computed masses in Table 1. 



Figure 13. Histograms of the cluster masses for the whole set 
of clusters (see Fig. 6) and of the sample of Table 2. 

—2.13 ±0.15 (r = 0.96) for our complete sample. These 
values compare with —2.2 given by Reddish (1978), who 
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included associations, too. The peak value of our mass 
function is at I26M0, instead of ~ 360 Mq of Reddish 
(1978). 

5. Conclusions 

In this paper we presented a standard method to syn- 
thesize in a reliable and automatic way the U,B,V lu- 
minosities of stars in an open cluster. One of the most 
important parts of the work is the selection of member 
stars based, as usual, on radial velocities and proper mo- 
tions (when available) and HR diagrams. We show how a 
convergence within 0.05 magnitudes of the integral mag- 
nitude is reached when all the stars in the HR diagram 
down to 6 mag fainter than the brightest member are 
included. 

The comparison of our photometric results for a sam- 
ple of clusters for which integral B and V data are found 
in Lynga's 5th Catalogue of Open Cluster Data indi- 
cates a high degree of compatibility between the two sets 
of data. This makes us confident about the use of the 
synthetic fluxes given by Lynga's catalogue (whose in- 
tegration procedure is not discussed) and, consequently, 
about our enlargement of that sample. An important new 
contribution is the evaluation of U magnitudes and U — B 
colours for 115 clusters. 

We also discussed a method to evaluate the luminous 
mass of a cluster via the computation of its cumulative 
mass distribution; this allowed us to calibrate a mass- 
luminosity relation which was used to obtain the present- 
day mass function of all the clusters with available integral 
My and age. The exponent of the high mass tail of the 
mass function is —2.04, and it represents an updating 
(being also based on a much more abundant sample) of 
Reddish's (1978) —2.2 result. This mass function is flatter 
than that of stars in the solar neighbourhood. A compar- 
ison of our mass estimates (referring to luminous stars 
in a cluster) with the dynamical estimates by Schmidt 
(1963) yields to a roughly constant ratio of non-luminous 
to the total mass in the cluster sample examined. Finally, 
the new integral magnitudes presented here lead to an 
increase up to 115 members of the set of clusters brighter 
than My = —4.5 and inside the 2 Kpc radius cylinder 
around the Sun, i.e. the number-complete sample defined 
in Battinelli & Capuzzo-Dolcetta (1991). This sample 
constitutes a natural luminosity -unbiased template for a 
comparison with star clusters in galaxies and gives precise 
constraints to evolutive models of the open cluster system 
in our galaxy. 
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